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The thermal quantum discord (QD) is studied in a two-qubit Heisenberg XXZ system with 
Dzyaloshinskii-Moriya (DM) interaction. We compare the thermal QD with thermal entanglement 
in this system and find remarkable differences between them. For instance, we show situations where 
QD decreases asymptotically to zero with temperature T while entanglement decreases to zero at 
the point of critical temperature, situations where QD decreases with certain tunable parameters 
such as D z and D x when entanglement increases. We find that the characteristic of QD is exotic 
in this system and this possibly offers a potential solution to enhance entanglement of a system. 
We also show that tunable parameter D x is more efficient than parameter D z in most regions for 
controlling the QD. 
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I. INTRODUCTION 

In recent years, the concept of entanglement [l| has re- 
ceived much attention in quantum information and quan- 
tum computation [2| due to its fascinating features in de- 
veloping the idea of quantum communication and quan- 
tum computer. 

Entanglement is one kind of quantum correlations [l[ 
and offers support for lots of quantum tasks such as 
quantum teleportation 0J3 i. quantum dense code [|| 
and quantum computing [a, lZ|- However, it is not the 
only part of quantum correlations So what is the 
quantum correlation without entanglement? An answer 
to this question is the notion of quantum discord (QD) 
given by Ollivier and Zurek 9]. It introduces a mea- 
sure of the quantumness of correlations (including en- 
tanglement). This notion has attracted much attention 
nowadays [j| IIol-[l5j . It can be treated as an important 
quantum resource and even used for the quantum com- 
putation task [T6rll8j . 

As a kind of quantum resource, entanglement can be 
generated, maintained and controlled in so many con- 
densed matter fields especially solid state systems such as 
quantum dots [H,|13|, Josephson junctions 21 1, trapped 
ions [22[ and cavity-QED systems (2^|, etc. Among sys- 
tems above, owing to its concise form of Hamiltonian 
and ability to be realized in real physical systems, there 
has been plenty of work discussing entanglement and 
thermal entanglement qualitatively and quantitatively in 
spin-chain systems [2J-129J. 

Since much work discuss the spin-spin interaction in 
spin-chain systems such as Heisenberg model, more and 
more researches study the influence of Dzyaloshinskii- 
Moriya (DM) interaction [3(1 E| which is arising from 
spin-orbit interaction in spin-chain systems [32H34j . Re- 



* Electronic address: 
t Electronic address: 



yxchen@zimp.zju.edu.cn 
zhiyin@zimp .zju.edu. en] 



cently, there is an interesting paper [35| investigating the 
thermal entanglement in an anisotropic XXZ model with 
DM interaction. The characteristic of thermal entangle- 
ment in that model is notable especially when you control 
the entanglement via tuning a;-component DM interac- 
tion. 

Since QD is considered as quantum correlation and 
quantum resource as well as entanglement. It is curious 
to ask what the characteristic of QD is at finite temper- 
ature, and what the differences are between thermal QD 
and entanglement in this system. In our paper we intro- 
duce this two-qubit anisotropic Heisenberg XXZ model 
with DM interation. We discuss the dependence of ther- 
mal QD in this system on temperature T and some other 
tunable parameters such as DM interaction parameters 
D z , D x . We also compare the QD with entanglement 
measured by concurrence [36H38j between two qubits in 
the model and find remarkable differences between them. 
QD is more robust than entanglement concurrence versus 
temperature T. (This property can also be found in Ref. 
[l4| where a Heisenberg XYZ model with a magnetic term 
was considered.) We find that the characteristic of QD 
is special in this system since QD decreases with certain 
tunable parameters such as D z and D x when entangle- 
ment increases. And this possibly offers a potential solu- 
tion to enhance entanglement of a system. Furthermore, 
we analyze the influence of different DM interaction pa- 
rameters D z and D x on QD and reveal that tunable pa- 
rameter D x is more efficient than parameter D z in most 
regions for controlling the QD. 

Our paper is organized as follows. In Sec. II we review 
the concept of quantum discord. In Sec. Ill we investi- 
gate the thermal quantum discord in two XXZ Heisen- 
berg model with different DM interaction parameters. 
We show the main results and compare thermal QD with 
thermal entanglement. In Sec. IV we discuss the opposite 
tendency of thermal QD and entanglement. Conclusion 
are then presented in Sec. V. 
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II. QUANTUM DISCORD 

In classical information theory, the correlation between 
two random variables A and B are described by mutual 
information which reads in two expressions. 

1(A : B) = H(A) + H(B) — H(A, B) (1) 

J(A : B) = H(A) - H(A\B) (2) 

Here, H(A) — — ^2 a Pa^og 2 p a is Shannon entropy. 
H(A, B) is the joint Shannon entropy which is defined 
by H(A,B) = -J2 a , b Pa,b^og 2 p a , b . H(A\B) is condi- 
tional entropy and introduced for quantifying the igno- 
rance (on average) about the value of A given B is known. 
This two expressions are equivalent using Bayes' rule: 
p(ai,bj) — p(ai\bj)p(bj) — p(bj\ai)p(a,i), where p{ai\bj) is 
the conditional probability. 

In order to generalize these expressions above to quan- 
tum case, we replace classical probability distributions by 
density matrices and the Shannon entropy by von Neu- 
mann entropy S(p) = — tr(plogp). Since we can denote 
the density matrix of a composite system AB by p AB 
and the density matrices of parts A and B by p A and 
p B , respectively, we can rewrite Eq. ([T]) and Eq. @ as 



S( P A ) + S( P B ) S( P AB ) 
J{ P A : p B ) = S(p A ) - S{p' 



~ A \P B ) 



(3) 
(4) 



where S(p AB ) is quantum joint entropy and S(p A \p B ) 
is quantum condition entropy. However, here we can not 
just obtain quantum condition entropy simply via replac- 
ing Shannon entropy by von Neumann entropy. Due to 
the definition of condition entropy, S(p A \p B ) depends 
on the choice of measurement. To get quantum con- 
dition entropy, we should choose a measurement which 
is described by a set of projectors {fif } performed lo- 
cally on B. After the measurement, the state of system 
is pi = fl B p AB U B / Pl where Pi = trftlf p AB fif ). Now 
we can get the definition of quantum condition entropy 
Following the Eq. ((¥]), another expression of quantum 
mutual information is 

J(p AB : {flf}) = S(p A ) - S(p AB \{fl B }) (5) 

Since different choices of {flf } decide different 3{p AB : 
{fif }), Eq. (]3]) and Eq. (|4|) may not be equal any more. 
According to the original definition Q , quantum discord 
is the minimum of difference between Eq. and Eq. 
©, i.e. 



D(p 



AB\ 



min Z(p A : p B ) - J(p AB : {nf}) 



(0) 



It is a measure of the quantumness of a pairwise corre- 
lation. With the expression of quantum mutual infor- 
mation of Eq. ([5]), following Ref. [9j, |39(, we can define 
classical correlation between A and B as 



CC{p AB ) 



max(J(p 
{nf} 



AB 



{nf})) 



(7) 



Eventually, we get quantum discord in another expres- 
sion: 



Q(p AB )=l(p A :p B )-CC(p AB ) 



(8) 



Here the quantum mutual information is used as a mea- 
sure of total correlations (iTJ, . 

Quantum discord is a measure to quantify all the quan- 
tum correlations including entanglement of a pairwise 
correlation. For two subsystems A and B are correlated 
classically only, the QD is zero. Moreover, for some pair- 
wise system such as Werner states @, QD is non-zero 
although the system is separable. In this sense, some two- 
body states that don't have entanglement do not mean 
there is no quantum correlations between them. We can 
even use this kind of quantum correlation power to do 
some quantum computation tasks [l6l - [l8j . 



III. XXZ HEISENBERG MODEL WITH 
DIFFERENT DM INTERACTION PARAMETERS 

In this section, we introduce an anisotropic XXZ 
Heisenberg model with different DM interaction param- 
eters. We show the dependence of thermal QD on tem- 
perature T and some other tunable parameters such as 
J z , D z and D x in this model. We compare the thermal 
QD with thermal entanglement concurrence in the same 
model and find some attractive results. 



A. XXZ Heisenberg model with DM interaction 
parameter D z 

The Hamiltonian H of a two-qubit anisotropic Heisen- 
berg XXZ chain with DM interaction parameter D z is 

H = Jo'ia* + Jo* a* + J t afa^ + D z (a%a% (9) 

where J and J z are the real coupling coefficients, D z is 
the z-component parameter of the DM interaction. Here 
the so-called DM interaction is a supplemented magnetic 
term arising from the spin-orbit coupling (30l . l3l| . a % 
[i = x, y, z) are Pauli matrices. In our paper, we consider 
the antiferromagnetic case (J > 0, J z > 0). This model 
is reduced to the isotropic XX model when J z = and 
to the isotropic XXX model when J z — J. Parameters 
J, J z and D z are dimensionless here. 

First we show the matrix form of this H in the standard 
basis {|00), |01), |10), 11)}, and get 



H = 



Jz 














— Jz 


2 J + 2iD z 








2 J - 2iD z 


-J z 














Jz 



(10) 



Now we want to consider the state of the system as 
thermal equilibrium state in a canonical ensemble. The 
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state of spin chain system at thermal equilibrium is given 



by the Gibb's density operator p(T) 



expj-pH) 



where 



Z = tr[exp(— 0H)] is the partition function of the sys- 



tem, H is the system Hamiltonian and f3 = 



K B T 



with 



T temperature. Kg is the Boltzmann's constant which 
we take equal to 1 for simplicity. Here p(T) presents a 
thermal state. Plenty of paper have discussed about en- 
tanglement in thermal state [241429} which is called ther- 
mal entanglement [42j. In our paper, we discuss another 
kind of quantum correlation-quantum discord at finite 
temperature i.e., thermal quantum discord. 

By some routine calculations, we can get expression of 
P{T): 



P(T) = - 



-pj z 







ve 





u 







ve l1 
u 










(11) 



where u = 

W = 



i(l + fAP^eW.-**) t 



1(1 



z 2pj z cos h(2/3u;)] and 



Z = 2e- pJ *[l -I 

D%,9 = arctan(^)- 
According to Eq. @ , to obtain QD of a two-qubit sys- 
tem, we should first get p A , p B for subsystems A and B 
respectively and p AB for the composite system. Since we 
have p AB already, we need to evaluate p A , p B . Choosing 
0), 1) as a basis, we can trace out the rest of the subsys- 
tem to get reduced density matrix of A and B employing 
formulas below. 



N B 

p A ^tr B p AB =£(/ai 
j'=1 



{^\)P AB {IA^ 3 )) (12) 



P B = tTAP 



N A 



iB) P AB m^i B ) (is) 



where I a and Ib are the identity operators, \<j>j) (j — 
1,2, ...,Nb) and \t/jj) (j = 1,2, ...,Na) are orthonormal 
bases in We and Ha respectively. 

Now we can get pa — Pb — j;(e~^ J * + u)I where 
/ is the identity operator. Next we evaluate quantum 
discord numerically according to the following instruc- 
tion. Choose the set of projectors {|V'i)(V'iU^'2)(^'2|}, 
where = cos6>|0) + e^sin6>|l) and \ip 2 ) = -cos0|l) + 
e _J * , sin0|O) , to measure one of the subsystems. The 
classical correlation CC (p AB ) is obtained numerically by 
varying the angles 8 and </> from to 2tt and eventually 
according to Eq. (JSj) we can get QD as well. 

Using the method in Ref. |36l - [38j . and results in [35j . 
we can also simply get the concurrence (In this paper, 
we use concurrence as the measure of entanglement) of 
this model at finite temperature. C{p{T)) = ^(e 2/3w - 
e -WJ* _ e -20w _ e -2/3J.) if j z > _ Wj and C(p(T)) = 

if J z < —w. 

Now by fixing some parameters we can see the vari- 
ations of thermal QD and entanglement under the in- 
fluence of other parameters. Looking at Fig. 1, where 



(a) 




(b) 




(c) 




FIG. 1: The QD (a), (c) and concurrence (b) versus T in the 
system with z-component parameter D z . Here J = 1. (a) 
Here J z = 0.2 and for the thick line D z — 1.0, normal line 
D z = 0.7, dashed line D z = 0.5. (b) Here J z = 0.2 and for 
the thick line D z = 1.0, normal line D z — 0.7, dashed line 
D z — 0.5. (c) Here J z = 1 and for the thick line D z — 1.0, 
normal line D z =0.7, dashed line D z =0.5. 



J = 1, we can find that the QD decreases with increas- 
ing of temperature T . This tendency is similar to the 
situation of entanglement as shown in panel (b) where 
the conditions are the same. Differently, we can see from 
panel (b) that thermal entanglement will vanish when T 
reaches certain point. And the larger parameter D z is the 
larger critical point is. The point denotes so-called criti- 
cal temperature of the system. This behavior is caused by 
mixing of the maximally entanglement state with other 
states. In other words, it means that with the increasing 
of temperature T, thermal entanglement of the system 
will be disappeared at critical temperature and there is 
no entanglement resource to use then. However, as we 
see in panels (a) and (c), thermal QD will not reach zero 
point. It acts asymptotically near zero when the tem- 
perature is high. And in the temperature regions where 
thermal entanglement is zero while QD is still alive. We 
even find that we adjust T very high (not shown in the 
plots), QD becomes very small but nonzero. This inter- 
esting phenomenon can be explained as follows. With 
the temperature grows, the role of thermal fluctuations 
exceeds quantum ones. However, it cannot kill the quan- 
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tumness. As we know, entanglement is not the only part 
of quantum correlations. Based on Ollivier and Zurck's 
argument || , it shows that absence of entanglement does 
not imply classcality. If you want to destroy the quan- 
tumness of a system, you should use the process of deco- 
herence. Decoherence will eventually make the quantum- 
ness of a system to be zero and then the system will be in 
totally classical state. Since QD measures total quantum 
correlations, it is obvious that even temperature is high, 
without introducing the decoherence, quantum correla- 
tions will not vanish. So in this sense, QD is robust than 
entanglement at a finite temperature. 

Focusing now at panel (a) of Fig. 1, when tempera- 
ture T increases, thermal fluctuation takes control of the 
system in this region, and the role of coupling constant 
such as J, J z and D z become weaker. So we can see the 
lines in panel (a) of Fig. 1 become concentrated together. 
Now we enhance the parameter J z to 1.0 as we can see 
in panel (c) of Fig. 1. We find that at the same temper- 
ature, lines become more separate than the case of panel 
(a). This effect is because of the influence of coupling 
constant J z makes the quantumness considerable. Next, 
let us take a look at panels (a) and (c) again and con- 
sider the influence of parameter D z to thermal QD. Due 
to the symmetry of DM interaction parameter, we can 
only take the strength of D z into account. So here only 
consider the D z > case. Apparently, we can see a re- 
markable phenomenon (Also shown in Fig. 3) in certain 
region that when the parameter D z increases, thermal 
QD decreases while this trend of QD is in contrast to the 
behavior of entanglement. Since QD and entanglement 
are both kinds of measure of quantum correlations. This 
is very wired. Worth to be mentioned, this phenomenon 
appears in D x case as well, and we will analyze this in 
detail in Sec. IV. 

Moreover, take a look at Fig. 2, where we fix J = 1 and 
D z — 1. Asymptotical phenomenon of QD also appears 
at finite temperature for different J z in panel (a). How- 
ever, we can see that QD increases with increasing of J z . 
This is corresponding to the case of entanglement shown 
in panel (b) . Here, J z as a anisotropic coupling constant 
of Heisenberg model, plays the same role in thermal QD 
as in the case of entanglement. 



B. XXZ Heisenberg model with DM interaction 
parameter D x 

Here we consider the case of the two-qubit anisotropic 
Heisenberg XXZ chain with DM interaction parameter 
D x . The Hamiltonian H' of this model reads 



H' = Jolol + Jo\ol + .J z o 7 io z 2 +D x (oXol~o\ol). (14) 

where D x is the x-component parameter of the DM in- 
teraction and other parameters are the same as in part 
A of this section. 



(a) 



0.30 
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0.15 
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FIG. 2: The QD (a) and concurrence (b) versus T in the 
system with z-component parameter D z . Here J = 1, D z — 1. 
(a) Here for the thick line J z = 3.0, normal line J z = 2.0 and 
dashed line J z = 1.0. (b) Here for the thick line J z = 3.0, 
normal line J z = 2.0 and dashed line J z = 1.0. 




FIG. 3: (Color online) The quantum discord is plotted versus 
T and D z where J = 1 and J z = 0.2. 



Following the exact same approach in part A, we can 



get 



H' = 



t J * 
-W x 

iD x 

V o 



iD x 









2J 


iD x 


2J 




-iD 


-iD x 


iD x 


Jz 



(15) 



^ = 2^ 




(16) 



where w' = y/(J + J z ) 2 + \D 2 X 
if = arctan( J+ 2 J ^ tu/ ), and 

M± = e -PJ, ± (g/3(J-u/) gin 2 ^ _, 



■ = arctan( J+J ^ w , ) 
MJ+ W >) gin 2 ^ 



■ + e 



£ = ie 



sm 



4> cos (f> + ie' 3( -' 7+t0 ' sin ip cos tp 



cos <p) 



Z' = 2e~P J cosh[/3( J - J z )] + 2e /3J cosh((3w'). 
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FIG. 4: The QD versus T in the system with z-component 
parameter D x . Here J = 1, .D^ = 1 and for the thick line 
J z — 0.9, normal line J z = 0.4, dashed line J z = 0. 



(a) 




(b) 




Also we can get p and p via straightforward calcu- 



lation: p' A = p' B = ^7 (a^ 
operator. 



v+)I where / is the identity 



Now, first we look at Fig. 4, where J = 1 and D x = 1. 
We can see that when J z increases, the thermal QD be- 
comes higher. This trend of QD is the same like entangle- 
ment case. (For better comparison we choose the exact 
the same value of parameters in Fig. 4 as in Ref. (35[) 
And also, this fits with the behavior of J z in the system 
with z-component parameter D z . 

While looking at Fig. 5, we find that besides the 
asymptotical phenomenon mentioned above, we can also 
find the opposite trend of thermal QD and entanglement. 
So here we can derive that different DM parameters can 
both influence the behavior of QD and entanglement and 
make them move towards different directions in certain 
region. We will address this in Sec. IV. 

In addition, back to the Fig. 1 and Fig. 5 which 
demonstrate the behaviors of thermal QD in aniosotropic 
Heisenberg model with parameters D z and D x . It is ob- 
vious there are some differences between two different 
cases. First, let us look at panel (a) of Fig. 1 and panel 
(b) of Fig. 5, we find that in the situation where other pa- 
rameters are exactly the same, D x is more reliable than 
D z when T becomes higher. This indicates that D z is 
easy to influence by the increasing of T (This also shown 
in Fig. 3 and Fig. 6). And from Fig. 1 and Fig. 5 we can 
find that if we tune the parameter J z larger, thermal QD 
depending on different D x change more obviously than 
the case of parameter D z . This means that parameter D x 
is more sensitive in certain condition along with change 
of parameters such as J z . So we can derive from this 
results that it is more efficient to use parameter D x to 
control the QD. In addition, we can also find that when 
tune parameter J z to 1.0, the two plots (panel (c) of Fig. 
1 and panel (a) of Fig. 5) become similar, which is mainly 
because the DM interactions now are not the dominating 
parameters controlling the behavior of the system. 



FIG. 5: The QD versus T in the system with x-component 
parameter D x . Here J = 1. (a) Here J z = 1.0 and for 
the thick line D z = 1.0, normal line D z = 0.7, dashed line 
D z = 0.5. (b) Here J z = 0.2 and for the thick line D z = 1.0, 
normal line D z =0.7, dashed line D z =0.5. 



0.0 T 
u " z 0.4 




FIG. 6: (Color online) The quantum discord is plotted versus 
T and D x where J = 1 and J z = 0.2. 



IV. DISCUSSION ABOUT OPPOSITE 
TENDENCY OF THERMAL QD AND 
ENTANGLEMENT 



In this section, we will discuss the opposite tendency 
of thermal QD and entanglement appeared in Sec. III. 

In previous section, we find that, in certain regions, 
when parameter D z increases, thermal QD decreases 
while at the same time thermal entanglement increases. 
This also appears in the case of parameter D x . So we 
say that when we tune the parameters of DM interaction, 
entanglement and QD have the opposite trend of move- 
ment. Obviously, this notable phenomenon is very wired 
since both QD and entanglement quantify the quantum- 
ness of correlations. There are some underlying reasons 
about that. Although QD and entanglement both kinds 
of measure of quantum correlations. Entanglement is a 
non-local correlation while QD figures a measure of whole 
quantumness of correlations for a pairwise system. There 
are some other quantum correlations without entangle- 
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ment which can be characterized by QD but not any 
kinds of entanglement measure such as concurrence. And 
we reveal that these other quantum correlations have re- 
markably different characteristic from the entanglement. 
Recent researches also suggest that these other 

quantum correlations can also be an important resource 
and maybe take respondence in mixed state quantum 
computation. Secondly, DM interaction is a magnetic 
term arising from the spin-orbit coupling. The staggered 
DM interactions bring on helical magnetic structures due 
to its form D • [S 1 x S 2 ], where Si, S 2 denote the spins 
and D is DM vector. DM interaction can contribute the 
strong planar quantum fluctuations which pose the align- 
ment ordering to both qubit-1 and qubit-2 in two-qubit 
system. Some researches [IH HiJ have shown that in cer- 
tain condition, DM interaction can build entanglement. 
So it is acceptable that entanglement increase with in- 
creasing of DM interaction in this paper. 

Now we can see that since QD figures the whole 
quantum correlations including entanglement, the phe- 
nomenon here indicate that, the other quantum correla- 
tions without entanglement we mentioned above decrease 
rapidly (more rapidly than the increase of entanglement 
at the same time) with increasing of parameters of DM 
interaction. So we can see from Fig. 1 and Fig. 5. 
that when QD decreases while the entanglement increases 
with the increasing of DM interactions. Here, this phe- 
nomenon can coincide what we have raised above that 
the characteristic of QD is exotic. Furthermore, using 
this characteristic, under certain conditions, we can tune 
the parameters of DM interactions to increase the entan- 



glement of a system and meanwhile decrease the other 
quantum correlations of the same system. The advantage 
of the method is obvious. If you want to take entangle- 
ment as a quantum resource to do quantum tasks, you 
want to enhance your entanglement resource and make 
the other quantum correlations of system as less as pos- 
sible, this approach would offer a possible solution. 



V. CONCLUSION 

The thermal QD of a two-qubit anisotropic Heisenberg 
XXZ chain with two distinct DM interaction parameters 
is investigated. We find that thermal QD is more ro- 
bust than thermal entanglement versus temperature T. 
And numerical results indicate that it is more efficient to 
control QD in the model with parameter D x than with 
parameter D z . Furthermore, with the increasing of pa- 
rameters D z and D Xl thermal QD varies oppositely to 
the case of entanglement. This phenomenon reveals that 
quantum correlations without entanglement have quite 
exotic characteristic from entanglement. In addition, this 
effect also offers a solution to make the entanglement of 
a system enhanced. 
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